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Abstract. Unified geometric approach to describing kinematics of elastic and plas- 
tic deformations of continuous media is suggested. On the base of this approach 
we study mechanical deformations, viscous flow, and heat transport in glassy plastic 
media. As a result we derive appropriate differential equations for these phenomena 
in a form applicable to liquids, elastic solids, and to plastic solid materials as well. 



1. Introduction. 

AFM-assisted^ electrostatic nanolithography (AFMEN) was suggested in [1] and 
[2]. It is new technique for nano-scale patterns formation on planar polymer films. 
By this paper we start series of papers aimed to prepare theoretical background for 
further numeric simulation of AFMEN process. First of all we consider kinematics 
(i. e. time-dependent geometry) of deformations and apply differential geometry to 
study it. Our main result is a method for separating elastic and plastic deformations 
within general nonlinear deformation tensor. Then we build this method into the 
standard framework of balance equations traditionally used to describe dynamics 
and thermodynamics of moving continuous media. Applying electrostatic field, 
surface phenomena, and phase transitions will be considered in separate papers. 

As far as our technique is concerned, we use curvilinear coordinates from the 
very beginning. This looks a little bit tricky, but in this way we make transparent 
tensorial nature of all quantities we use. Moreover, the problem of AFMEN pro- 
cess simulation, which we are going to study numerically, has obvious cylindrical 
symmetry. So we prepare use of cylindrical coordinates in our future calculations. 
For reader's convenience in section 2 just below we resume some well-known facts 
concerning curvilinear coordinates. 

2. Moving frame of curvilinear coordinates. 

Let x^, x^, be three Cartesian coordinates of a point. Below we use both up- 
per and lower indices, following traditions of tensorial analysis (see [3]). Curvilinear 
coordinates j/^, y^, of a point are usually introduced by functions 

x^ = x^{y^,y'^,y^), 

x^ =xHy\y^y^), (2.1) 
= x^{y^,y'^,y^). 



^Here AFM is abbreviation for Atomic Force Microscope. 
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Functions (2.1) define transition to curvilinear coordinates. Transition back to 
Cartesian coordinates is determined by similar three functions 

< x^=xHy\y^,y^), (2.2) 
^ x^ =x^{y\y'^,y^). 

Maps (2.1) and (2.2) are inverse to each other. Therefore if we define their Jacobi 
matrices 5* and T by partial derivatives 



Si = 1^, Tl = 1^, (2.3) 

they are also inverse to each other: T = S~^. Matrix S is called direct transition 
matrix, while T is called inverse transition matrix. 

Let ei, 62, 63 be constant frame of Cartesian coordinates ^, i.e. these 

arc three base vectors directed along three coordinate axes. Then, using functions 
(2.1), we can define vectorial function r = r(j/^, y^, y^): 

3 

v^Y.^'{y\y\y^)ei. (2.4) 
1=1 

This is radius-vcc;tor of a point expressed through c;TU'vilinc!ar coordinates. Differ- 
entiating (2.4) with respect to y^, y^, y^, we get three vector- functions: 

Ei = Ei(yi,y2,y3) = 0, i = 1,2,3. (2.5) 

Vectors Ei, E2, E3 form a frame of curvilinear coordinates y^, y^, y'^. This is mov- 
ing frame, since vectors Ei, E2, E3 depend on coordinates of a point to which they 
are attached. Transition matrices (2.3) relate moving frame Ei, E2, E3 to constant 
frame 61, 62, 63 of Cartesian coordinates and vice versa: 



E, = f^5/e„ 6,=X:r/E,. 

g=l q=l 

Mutual scalar products of frame vectors (2.5) define fundamental tensor of our 
space. This is metric tensor g with components: 

Qij = (Ei, E,). (2.6) 

Matrix (2.6) also is known as Gram matrix of moving frame Ei, E2, E3. Inverse 
matrix to (2.6) define another fundamental tensor. This is dual metric tensor, by 
tradition it is denoted by the same letter g, but for its components upper indices 
are used: g'^K Being mutually inverse, metric tensors are related by formula 



9=1 
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Here are Kronecker symbols. They determine components of unit matrix: 



1 for i = j, 
for i ^ j. 



3. Tensor fields in curvilinear coordinates. 

All tensor fields in curvilinear coordinates y^, y^, y^ are referenced to moving 
frame Ei, E2, E3. This determine some features of their differentiation. Thus, if 
(fi = (f{y^,y'^,y^) is a scalar field, then applying gradient operator V to it we get 
three components defined by partial derivatives: 

Vi^=|^, i = 1,2,3. (3.1) 
dy' 

However, unlike (3.1), applying gradient operator to tensor field X is more compli- 
cated procedure. For components of resulting tensor field VX we have formula 

■ - dx'^-Y 

^^^h-js - Qym 

r 3 s 3 (3-2) 
I F^fc vn--- Ofc--- «r _ pf'fc Y*i 

fe=lafe = l fe=l 6j,=l 

Here T^j is three-dimensional array of connection components or Christoffel symbols. 
They are defined by metric tensor g according to the following formula: 

■' 2 \ oy^ oy^ oy^ ' 



Formulas (3.2) and (3.3) arc well-known in differential geometry (see [3]). 

In addition to metric tensor (2.6) and connection components (3.3) there are 
also volume tensor u: and dual volume tensor denoted by the same symbol: 



LOijk = \/detg £yfe, w'-''' = —7==, (3.4) 

Vdetg 

Here Sijk = ^'■''^ are Levi-Civita symbols. They are defined as follows: 



if i = j, j = k, OT k — i; 

1 if {ijk) is even transposition of (123); (3.5) 
— 1 if {ijk) is odd transposition of (123). 



Levi-Civita symbols (3.5) do not form a tensor. However, supplying scalar factors 
to them, we get two tensor fields (3.4). 

4. Deformation of continuous medium. 



Suppose that our space is filled with medium of some kind. Deformation of 
medium is due to the displacement of its points. Suppose that the point with 
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coordinates y^, y^, has moved to the point with coordinates y^, y^, y^. This 
situation is expressed by the following three functions: 



2 - y'^{t,y\f,f), 



y 



(4.1) 



^ ~y\t,y^,f,y^). 



Argument t in (4.1) is responsible for time evolution of displacement. Time deriva- 
tives of these functions determine velocity vector components: 



v'=y' = ^, i = 1,2,3. 



(4.2) 



Velocity vector v itself is calculated as a sum representing its expansion in moving 
frame at the point with coordinates y^, y^, y^: 



(4.3) 



Functions (4.1) define time-dependent map from space to space. Let's denote it by 
r. Then inverse map is given by similar functions 



y^ = y^{t,y^,y'^^y^), 
y^ = y'^{t,y\y'^,y^), 
. y^ = y^{t,y\y'^,y^)- 



(4.4) 



Quantities v', as defined in (4.2), are functions of coordinates y^, y^, marking 
initial position of a point of medium. Using (4.4) we can convert them to the 
functions of coordinates y^, y"^, y^ marking current actual position of that point: 



Time dependent maps (4.1) and (4.4) define two Jacobi matrices S and T: 



(4.5) 



qi _ 9y' 
^ dyi 



fi _ 9y' 
^ dyi ' 



(4.6) 



Matrices (4.6) are inverse to each other and quite similar to that of (2.3). However, 
in contrast to matrices S and T, these two matrices depend on t and describe 
physical state of our continuous medium (more precisely, they describe the state of 
deformation). Matrix T has proper arguments t, y^. y^, y^, while arguments of S 
should be corrected by inverse map (4.4), i. e. S ^ S ° t~^. 

Let's take vectors of moving frame Ei, E2, E3 at the point with coordinates 
y^, y^, y"^ and send them to the point with coordinates y^, y^, y^ by means of map 
T~^. As result we get another frame Ei, E2, E3: 



1,2,3. 



(4.7) 



r=l 



NOTE ON KINEMATICS, DYNAMICS, AND THERMODYNAMICS ... 5 
Mutual scalar products of frame vectors (4.7) form a matrix with components 



3 3 



Gi^ = {%, E,) = ^ ^ 9rs{y\y\f) f[ f/. (4.8) 

r=l s=l 

Upon transforming all arguments in (4.8) to t, y^, , we get tensor field G 
with components Gij = Gij{t, y^, y^, y^). Tensor G at the point with coordinates 
y^, y"^, y^ is an exact quantitative measure of deformation of our medium at that 
point. For small deformations described in Cartesian coordinates we have 

Gij{t, x^,x^,x^)= Qij -2 Uij + ... (4.9) 

Tensor u in (4.9) is standard tensor of deformation as defined in [4]: 

1 / dui duj \ / , , „x 



2 \dxi dx\ 

By dots in (4.9) we denote terms of higher order with respect to small displacements 
= fe^, V? = Sx^, = dx^. Relying upon above formulas (4.9) and (4.10), now 
we define deformation tensor u as follows: 

Uij = ^^4^. (4.11) 

Using (4.11) instead of (4.10), we have not to restrict ourself to small displacements 

and can consider deformations of any magnitude in any curvilinear coordinates. 

Let's calculate time derivative for deformation tensor u defined by formula (4.11). 
Differentiating (4.8), by direct calculations we derive the following formula: 



Gij = - E VfeGi,- v''-J2 ^iv" - V,•^;^ (4.12) 

k=l fe=l fe=l 

Now, combining formulas (4.11) and (4.12), we obtain ultimate expression for inf. 

3 

Uij = - ^ {VkUij v'^ + u^j Viv"^ + Vjv'^) . (4.13) 



Let's denote by Vij first term in right hand side of (4.13): 

V, = ^n^. (4.14) 



Last three terms under summation in formula (4.13) are nonlinear with respect 
to deformation functions (4.1). They are usually omitted in the case of small 
deformations. In that case we would have iiij = Vij. 



6 



S. F. LYUKSYUTOV, R. A. SHARIPOV 



5. Dynamics of continuous medium. 



Dynamics of any continuous medium (either liquid, solid, or gaseous) is usually 
described in terms of three balance equations. These are equations for 

(1) mass balance; 

(2) momentum balance; 

(3) energy balance. 

Mass balance is most simple among balance equations. It is written on the base 
of the following statement; change of mass enclosed in any fixed volume within 
continuous medium is determined by mass flow through its boundary: 



Sum in (5.1) is divergency divj, where j = pv is density vector for mass flow. 
Components of velocity vector are determined by formula (4.5), while operator Vfc 
in (5.1) should be applied according to the formula (3.2). 

Momentum balance equation is more complicated, though it has the structure 
similar to mass balance equation (5.1). Momentum balance is written as 



Vector f with components p, p in right hand side of (5.2) determines density 
of volume forces in continuous medium. Symmetric tensor 11 with components n**^ 
determines density of momentum flow in continuous medium. Exact formulas for 
tensors f and 11 depend on various properties of medium and on those phenomena 
we are going to consider. 

Energy per unit of volume in continuous medium is a sum of two components: ki- 
netic energy and inner thermal energy due to chaotic motion of atoms and molecules. 
Therefore energy balance equation is written as follows: 



Vector w with components w^, w^, in (5.3) determines density of energy flow. 
Scalar e determines energy production/dissipation per unit volume of medium. It 
is due to work of force f in (5.2) and from other possible sources (e. g. Joule heating 
due to electric current in conducting medium). Potential energy is not included 
into left hand side of (5.3). However, work of potential forces is taken into account 
among other terms as a part of scalar e in right hand side of (5.3). 

Deformation state of a medium is completely determined by tensor G and metric 
tensor g. Prom (4.12) and (5.1) one easily derives: 




(5.1) 




(5.2) 




(5.3) 



Inp — 



In det G — In det g 



= const . 



(5.4) 
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Tensor 11 in (5.2) is usually given by the following formula: 

n''' = pv'v'' - a'''. (5.5) 
Here a'^ is stress tensor. Substituting (5.5) into (5.2), we get 

+ ^ V. ( p .*.'=)= r + ^ V^a^^ (5.6) 
fc=i fe=i 

Now, taking into account (5.1), from (5.6) we derive 

k=l k=l 

Using (5.4), now we can calculate time derivative for the density of kinetic energy: 

l(^) + E V. (^^ .') - 1 /* + E E «. v..^" 

^ ^ fe=l ^ ^ i=l i=l fe=l 

Looking at formula (5.8), we see that right hand side of this formula is completely 
determined by parameters ct*'^ and /* from (5.5) and (5.2). 

6. Elastic solids state media. 

Deformation state of a medium is completely determined by tensor G and metric 
tensor g. Therefore e in (5.3) satisfies the following equality: 

de = Tds-±±'-^. (6.1) 

Here s is an entropy per unit mass of solid state medium. Note that ct'^ in (6.1) 
should not coincide with a^^ considered above. Free energy per unit mass is deter- 
mined by standard formula f = e — T s. Therefore 

df = -sdT-±±'-^. (6.2) 
i=i j=i ^ 

Free energy per unit mass / is a function of temperature T and deformation state 
of medium: / = /(T, G). From (6.2) we derive 



2p dG, 



(6.3) 
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Solid materials can exhibit different properties in different directions. We consider 
only those polymer materials, which are homogeneous and isotropic (i. e. uniform 
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in all directions). For such materials /(T, G) is given by formula 

/=/(T,A[i],A[2],A[3]), (6.4) 
where A[i], A[2], A[3] are three scalar invariants for linear operator G: 

_tr(G) _ tr(G-G) _ tr(G • G • G) 
'^[i] - — ^> \2] ^ , A[3] . (6.5) 

Linear operator G in (6.5) is determined by its matrix Gj, where 

3 

G;=J29''Gkj. (6.6) 

k=l 

Using (6.5) and (6.6), we calculate partial derivatives 

dG-r^' 8G-r^^ dG^^^hu 

Applying (6.3) and (6.7) to (6.4), we find most general formula for a^^: 

3 3 

= /[I] g'' + fm G''^ fm g'" gu, , (6.8) 

fe=l g=l 

Here /[ij, /[2], /[sj are coefficients depending on T and on scalar invariants (6.5): 
f^^=-^^^J^^2h^lMM, .^^^2,3. (6.9) 

Formula (6.8) is exact, but very huge. Its use in numeric simulation is impossible 
not because of huge computations, but since function (6.4) is not properly measured 
experimentally for broad range of its arguments. Formulas (6.8) and (6.9) are worth 
for us only because they indicate the dependence 

a'^ =a'^{T,G). (6.10) 

Traditionally linearized version of this dependence (6.10) is used for the case of 
small deformations when Uij are much less than gij. 

Now let's consider formula (6.2) again. If function (6.4) is known, entropy per 
unit mass can be calculated as partial derivative: 

s = -%- (6.11) 

Duo to (6.11) we can treat s as a function s — s{T, G). Moreover, s depends on G 
through scalar invariants (6.5) just like function / itself: 

s = s(T,A[i],A[2],A[3]). (6.12) 
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Then for thermal energy per unit mass we get 

s=f + Ts = e(T, G) = £(T, A[i], Ap], A^]). (6.13) 

If parameters w;'^ and e in (5.3) are known, then, substituting (6.13) into (5.3), we 
get the equation for temperature function T = T{t,y^ ,y'^ ,y'^). 

In thermodynamics specific thermal energy e is often treated as a function of 
entropy. Indeed, inverting the dependence of s on T for fixed A[i], Apj, Apj in 
(6.12) and substituting T = T{s, A[i], A[2], A^j) into (6.13), we get 

£ = e{s,G) = £(s,A[i],A[2],A[3]). (6.14) 

In the absence of heat transfer and viscosity, parameters w'^ and e can be derived 

from entropy balance equation. In this case dynamics of solid state medium is 
adiabatic. Therefore one can write the equation 

^ + ^V.(p..'^)=0. (6.15) 

By analogy with (6.15), now we calculate the following quantity: 

^ + V. ( p £ .*=) = p I + p v.£. (6.16) 
fe=i fe=i 

Applying (6.1) and (6.14) and using Vkgij = 0, which is basic property of metric 
tensor, we get formulas for de/dt and Vfc£ in right hand side of (6.16): 

de _ ds -A a'^ Gjj 
dt~ dt 2p ' 

V.£ = TV..-EE^^- 
Substituting (6.17) into (6.16) and taking into account (6.15) and (5.1), we get 

^ + V.( p£ .'^) = - E E ^ fc, + ^ .'^ V,G^ . (6.18) 

fe=l i=l j=l \ fe=l / 

Now let's substitute (4.12) into (6.18). As a result of simple calculations we find 



dips) 



3 3 



^Vfc(p£^'=) = EE^^ Ev,«'=G,, + ^V,z;'=G 



fe=l i=l j=l \ fe=l fe=l 

3 3 3 /^J —fez ^3 3 3 / ^ 

= E E E ' = E E V.-. E 

i=l j=l fe=l i=l fe=l \ j=l y 
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Note that the above formula is quite similar to formula (5.8). Adding these two 
formulas, we obtain the following equality: 



■I l-'-^+pej +2^Vfc( ^-^v'' + pev'' 



dt\ 2 ' ) ' ^ 2 

3 3 3 ' 3 3 / 3 \ ^'^■^^^ 

1=1 j=i \j=i / 

Right hand side of the equality (6.19) simplifies substantially if stress tensor a^'^ 
and tensor ct*-' introduced in (6.1) are related as follows: 

3 

aik^Y^Q^-3k_ (6.20) 
Due to the relationship (6.20) formula (6.19) transforms to the following one: 

|2 \ 3 /„l,,|2 



= E^ir+EEVfe(^;i<7''=). 



3 3 (6.21) 



i=l i=l k=l 

The relationship (6.20) is self consistent. Indeed, if (6.20) is fulfilled then we have 
the equality (6.21) with quite transparent interpretation. First two terms in (6.21) 
describe energy increment and energy flow due to mass transport. Two terms in 
right hand side of (6.21) describe energy creation due to external forces and due to 
stress forces in continuous medium. Comparing (6.21) and (5.3), we get 

I 1 2 3 3 

w^ = ^^v''+pev''-Y,Via'\ e = Y,Vif- (6.22) 

i=l i=l 

Due to special form of tensor a and due to (6.20) stress tensor <t is symmetric: 
a^i = a^^. Indeed, from formula (6.8) for components of stress tensor we derive 

3 3 



a^i = /[ij G^' + E E 9k, 



'='^=' (6.23) 

3 3 3 3 ^ ' 

+ miY.T. /[31 9k, 9mn . 

k=l g=l m=l n=l 

Remember that coefiicients /[ij, /[2], /[a] in (6.23) are determined by formula (6.9). 
They depend on deformation G and temperature T. 

7. Heat transfer, viscosity, and entropy production. 

Formulas for 11*'^ and w'^ in balance equations (5.2) and (5.3) become more 
complicated if we take into account viscosity and thermal conductivity of medium. 
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In the case when Vv 7^ different parts of medium immediately adjacent to each 
other move with different velocities. This gives rise to forces of viscous friction. 
These forces are described by additional term in (5.5): 

W'' = pv'v'' -a''' - a'''. (7.1) 

Here cr'*^ are components of viscous stress tensor a. In linear approximation they 
are linear with respect to velocity gradients: 

3 3 

^^'^EE^^'''^^^- (7-2) 

j = l g=l 

Here Vjq are components of symmetric tensor introduced in (4.14), while ry**^-'? in 
(7.2) are components of viscosity tensor. They possess the following symmetry: 

^ikjq _ ^jqik _ ^kijq _ jjikqj _ (7-3) 

Components of viscosity tensor (7.3) are kinetic coefficients. In near equilibrium 
deformations they are functions of temperature and deformation tensor: 

^ifejg =^ifej9(r^G). (7.4) 

If continuous medium is non-uniformly heated, i.e. VT ^ 0, this may cause 
direct heat transfer without mass transport in it. This phenomenon is due to 
thermal conductivity of medium. Thermal conductivity of medium is described by 
additional term in formula for density of energy flow w: 

I |2 ^ 3 3 

^k ^ ^fc + p £ _ ^ ^. ^^fc _ ^ ^. _ y^y ^ik _ (75^ 

i=l i=l i=l 

As compared to (6.22), in (7.5) we have two extra terms. First is due to viscous 
stress tensor cr, second term contains components of heat conductivity tensor x**^. 
Like ri^^i'^ in (7.4), components of heat conductivity tensor are kinetic coefficients, 
they depend on temperature T and deformation G: 

x*'^ =x''=(T,G). (7.6) 

Heat conductivity tensor (7.6) is symmetric, i.e. x*'^ = >f*^'. 

Let's recalculate entropy balance equation (6.15), taking into account additional 
terms in (7.1) and (7.5). Instead of (6.21) now we have 

|(£M!,,.),gv.(^.>..„.). 



= E + E E '''' + ^'') + E E ^'^(^^^ ^'')- 



i=l i=l fe=l i=l fe=l 
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Equations (5.7) and (5.8) are replaced by the following two equalities respectively: 

J! 1 ' J^„1 ' 7^ 1 



fe=l ^ fe=l k=l ^ 



(7.8) 

3 3 3 3 3 ^ ' 

= E + E E + E E ^» ^fe^'"- 

i=l i=l fe=l i=l fe=l 

Subtracting (7.8) from (7.7) and taking into account (5.1), we obtain 

§ + E -^^^ ^'^^ = E E (V'^^' (^*' + ^'') + ^'^(^^^ ^'')) • c^-Q) 

fe=l i=l fe=l 

derivatives de/dt and Va;£ in left hand side of (7.9) we can apply (6.17). Then 

3 3 3 _jj / 3 \ 

fe=i j=i 3=1 \ fe=i / (7 \ n\ 

3 3 • J 

= EE (^-^^^ (a'^ +a''=) + Vfe(V.r^^*'=)) . 

i=l fe=l 

Now let's substitute (4.12) into (7.10) and take into account formula (6.20): 

+ E v^(p^^^) = EE +^^(^-^-") = 



P 

For 



dt 



k=l i=\ fc=l 



i=l k=l ^ ^ i=l k=l j=l q=l i=l k=l 



In ultimate form this is entropy balance equation generalizing equation (6.15): 

dips) A„ / ^ v,rx''=\ 

fc=l \ i=l / 

1=1 fe=l j=l q=l j=l fe=l 

Last two terms in (7.11) are positive. This means that total entropy is ever growing, 
when medium is evolving toward thermodynamic equilibrium. 

8. Liquid state media. 

Balance equations in liquids are quite similar to those in elastic solid materials. 
The only difference is that liquids are always isotropic. Deformation state of liquid 
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material is completely determined by its density p. Therefore, instead of (6.13), for 
specific thermal energy per unit mass we have 

e = e{T,p). (8.1) 
Like £ in (8.1), specific free energy per unit mass is given by formula 

f = f{T,p). (8.2) 
Density p is related to deformation G by means of formula (5.4). Hence 

dp d(lndetG) tr(G-i dG) ^J^G'UGij ^. 

J = 2 = 2 = 2.2.^—' 

where G^^ are components of inverse matrix G"^. Differentiating (8.2) and taking 
into account (8.3), due to the equality (6.2) we have 

=-p^^&^. (8.4) 
op 

Substituting (8.4) into (6.20), we obtain formula for stress tensor in liquid medium: 

a*^=V|^5'^- (8.5) 
Scalar factor in (8.5) is interpreted as pressure. Indeed, we have 

P = p'j-p^ ^''=-P9''- (8.6) 

Viscosity tensor in liquids simplifies and takes the following form: 

nil 



(8.7) 



Heat conductivity tensor in liquid medium is also simpler than in solid media: 

^ik=^gik_ (8.8) 

Pressure p in (8.6) and scalar parameters rj, and >r in (8.7) and (8.8) all are 
functions of temperature T and density p of continuous medium: 

P = P{T,P), r] = r]{T,p), 

C = aT,p), K = x{T,p). 



Other parameters in balance equations (5.1), (5.2), (5.3) for liquids are same as for 
solid state media. 
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9. Plastic materials. 

Saying plastic materials we mean pitch-like very dense sticky liquids and many 
solid materials with no crystalline grid, e.g. glass and polymer materials. They 
resist to deformations like solids and can flow like liquids, though sometimes very 
slowly. In order to describe such materials mathematically we need to divide de- 
formation G into two parts: elastic deformation G and plastic deformation G: 

G'i. = EEG,'=G,,G/. (9.1) 

fe=l 9=1 

Like Gij, elastic deformation tensor Gik in (9.1) is symmetric: Gik = Gki- Plastic 
deformation tensor Gkj is not necessarily symmetric. 

Plastic deformation arises as a response to stress tending to relax this stress. 
Elastic deformation tensor G is thermodynamic parameter of plastic medium, while 
plastic deformation tensor G is kinetic parameter. Therefore all thermodynamic 
quantities and kinetic coefficients for near thermodynamic equilibrium deformations 
of plastic medium depend on elastic deformation tensor G and on temperature: 

e = e{T,G), f = f{T,G), 

^ikjq ^ ^ikjq G) , h'^ = h'^ (T, G) . 

Components of total deformation tensor Gij in (9.1) satisfy differential equations 
(4.12). These equations can be written as follows: 

Q„ 3 3 3 

^ + E V.Gi,- = - ^ WiV' Gri - E Gir "^jV". (9.3) 
r=l r=l r=l 

For components of elastic deformation tensor Gij we write analogous equation: 

i~j 3 3 3 

_|i + ^ v.Gfc, = - ^ Vfci;'- Grq - E Gkr V,t;'- + 2 6^,. (9.4) 

r=l r— 1 r— 1 

Components of symmetric tensor in (9.4) are kinetic coefficients. Therefore 

efcg = efcg(T,G). (9.5) 

Now we need some special facts concerning all three deformation tensors in (9.1). 
Symmetric matrix Gij determined by formula (4.8) is non-degenerate and positive. 
Due to (9.1) we have the equality det G = det G • (det G)^. Hence 

det G 7^ 0, det G 9^ 0. (9.6) 



Moreover, due to (9.1) and (9.6) matrix Gkq is also positive. Now let's use formula 
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Gi=Y,9"'Gkj. (9.7) 



fe=i 



Theorem 9.1. For any symmetric matrix Qkq o,nd for linear operator G deter- 
mined by formula (9.7) there is unique symmetric matrix 6kq such that 



'^Qkq = Y.ekrGl + Y,Gl6rq. (9.8) 

r=l r=l 

This is purely mathematical fact with rather simple proof. Applying theorem 9.1 
to matrix (9.5) we get tensor 6 with components 

ekq = ekq{T,G). (9.9) 

Raising index in (9.9) we get symmetric linear operator 6 with components 

3 

e^^Y^a^'Ou,. (9.10) 
fc=i 

Due to (9.8), (9.9), and (9.10) we can write equation (9.4) in the following form: 

(9.11) 



Gkr Vgt'^ + O^. Grq + Gkr 0^ . 



r=l r=l r=l 

Now we are able to write differential equations for plastic deformation tensor: 

fif^k 3 3 3 

^ + ^ z;'- WrG^ = i^i - ^^'"^ ^r) - E Gl (9.12) 

r=l r=l r=l 

Forgetting principle. This is basic principle characterizing plastic deforma- 
tions that wc c;onsidcr in present paper. Suppose that plastic medium evolves from 
initial state with no deformation at time instant t = to some intermediate state 
Bitt = to, then it continues its evolution for t > to. Forgetting principle states that 
if in intermediate state total deformation of medium is purely plastic, then further 
evolution of medium will be so as if in intermediate state it had no deformation at 
all. Equation (9.12) is written on the base of this principle. It is compatible with 
(9.1) and with equations (9.3) and (9.11). 
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10. Thermodynamics of plastic medium. 

Balance equations (5.1), (5.2), and (5.3) remain unchanged for plastic medium. 
We also keep unchanged formulas (7.1), (7.2), (7.5), and second formula in (6.22). 
However, due to (9.2) and (9.11) we should revise formulas (6.2), (6.3), (6.4), (6.5), 
and (6.6). Formulas (6.2) and (6.3) are replaced by the following ones: 

For homogeneous and isotropic plastic material function f(T, G) is more special: 

/ = /(r,A[i],A[2],A[3]). (10.2) 

Though A[i], A[2], A[3] in (10.2) are different from that of (6.4), they are scalar 
invariants of linear operator G that was defined above by formula (9.7): 

tr(G) tr(G-G) tr(G • G • G) 

'^[l] = -^[2] = ^ ' ^[3] = ^ • (10-^) 

Next two formulas are analogous to (6.8) and (6.9): 

3 3 

= /[i] 9^' + /[2] G^' + E E /[31 9kg G"' , (10.4) 

fe=l g=l 

f - 2ip i9/(T,A;i],A[2],A[3]) -190 nn>^^ 
^[^1-"^ dX^^ ' *-l'2'3. (10.5) 

As we said above, we keep unchanged energy balance equation (5.3) with w'^ 
given by formula (7.5) and e given by second formula (6.22). This means that we 
keep unchanged formula (7.7). However, thermal energy e now depends only on 
elastic part of deformation tensor. Therefore, we should recalculate relation of a'-' 
and stress tensor cr*-'. Due to (10.1) and consequent formulas (10.2), (10.3), (10.4), 
and (10.5) instead of (6.14) now we should write 

£ = £(S,G) = £(s,A[i],A[2],A[3]). (10.6) 

By the same reason instead of formulas (6.17) we should write 
de „ ds a^^ dGij 



dt dt ^^^2p dt 



=1 j= 

3 3 



(10.7) 



Vfc£ = TVfcS-EE 



Applying (5.1) and (7.8) to (7.7), we derive the following equation: 

Q 3 3 3 

P-^+IZp"''^''^ = EE (^^^^^ (^'' + a^') + Vfe(V,T .<*'=)) . (10.8) 

fe=l i=l fe=l 
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Note that (10.8) is just the same as the equation (7.9). Now we should substitute 
(10.7) into this equation. As a result we obtain 

3 3 ^ ' 

i=l fe=l 

Applying (9.4) to the equation (10.9), we transform it to the following one: 
ds „ v^AAct'^' 



fe=i 1=1 j=i fe=i 

3 3 



i=l fe=l 

Now let's recall formula (6.20) and write analogous formula in present case 

3 

(,ik=^Q^^Jk_ (;^0.11) 

Applying (10.11) to (10.10), we simplify it substantially. Here is resulting equation 

^^+E^^^(/'-') = EE-^^%+ 

'='^=' (10.12) 
+ E E (^'^^^ ('^'^ + ^'') + Vfc(ViT x''^)) . 

i=l fc=l 

For &ij in (10.10) wc used formula (9.8). By further transformations similar to 
that of section 7 we can bring the equation (10.12) to the form analogous to (7.11): 

^.Ev.U-E^)^EE^. 

^ ' (10.13) 

' E E E E ''''' ^'''^^ I E E """^ ^^'^ 

As compared to (7.11) in (10.13) we have one extra term in right hand side. It 
should be positive hke other two terms in right hand side of (10.13): 

EE^^o- (10-14) 
i=l j=l 

Formulas (10.13) and (10.14) yield physical interpretation of tensor (9.9). This 
tensor determines entropy production due to plastic deformation of medium. 
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Formula (10.11) relates tensor cr and stress tensor a. Applying (10.11) to (10.4), 
we derive formula for stress tensor of isotropic plastic medium: 



3 3 



^io = /j^j + ^ ^ G^>^ ff,, G« 



(10.15) 

3 3 3 3 ^ ' 

+ E E E E /[3] G«™ gmn G"^' . 

fc=l g=l m=l n=l 

Formula (10.15) is analogous to (6.23). Like formula (6.23), this formula means 
that stress tensor <j is symmetric: cr*-' ~ a^^. In general case for non-isotropic 
media symmetry of tensor 11 in (7.1) and hence symmetry of cr is derived from 
conservation law for angular momentum (see [4]). 

Note that in (9.2) we declared that tensors cr and & depend only on elastic 
part of deformation tensor (9.1). However, for coefficients /[ij, /[2], /[3] in formulas 
(10.4) and (10.15) we have explicit expression (10.5) containing entry of density p. 
Due to (5.4) density of medium depends on total deformation tensor G, but not 
on its elastic part G only. In order to avoid this discrepancy we need to introduce 
additional restriction for plastic part of deformation tensor 

detG=l. (10.16) 

Then due to (9.1) and (10.16) we get detG = detG and (5.4) is replaced by 
analogous relationship binding density p to elastic deformation tensor G: 

In det G — In det g . ^ , „n 

Inp ^=const. (10.17) 

Restriction (10.16) and the equality (10.17) following from it are reasonable from 
physical point of view. Indeed, according to forgetting principle (see above) , plastic 
deformations are those which can be forgotten. However deformations changing 
density of medium cannot be forgotten. Hence they cannot be purely plastic. 

Restriction (10.16) leads to the restriction for tensor 6 in (9.12). Indeed, differ- 
entiating (10.16) and applying equation (9.12), we derive 

3 

tie = "^6^ = 0. (10.18) 

k=l 

Thus, formula (10.18) means that tensor determines symmetric linear operator 
with zero trace. It is symmetric due to (9.10) and symmetry 0i.q = Oq^. 

Plastic deformation is a way for draining stress. Tensor determines the rate of 
such draining. In elastic solid materials total deformation is purely elastic. There- 
fore Gj = 6j. Substituting = Sj into the equations (9.12), we find 

ej = 0. (10.19) 

Thus, elastic solid medium can be treated as limiting case of plastic medium with 
vanishing tensor 0? — > 0. 
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In liquids tensor is undetermined. Indeed, above in section 8 we treated liquid 
state media in the same way as purely elastic solid state media, but with special 
form of free energy function (see formula (8.2)). Then 9 is equal to zero like in 
(10.19). However, we could treat liquids as plastic media by introducing some 
arbitrary tensor G. In any case stress tensor in liquids is determined by formula 

a'i = -pg'i, (10.20) 

see (8.6) above. Substituting (10.20) into left hand side of (10.14), then taking into 
account (9.10) and trace condition (10.18), we derive 

i—l j—1 i—1 J — 1 i=l 

This means that tensor 6 for liquid media is not thcrmodynamically fixed. There- 
fore asymptotical behavior of functions O^g = 9kq{T, G) near phase transition point 
from plastic solid state to liquid state T — > Tp^.tr. requires separate investigation. 
This will be done in separate paper. 
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